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O " ABSTRACT 

z 

We study the formation, evolution and physical properties of accretion disks 
formed via wind capture in binary systems. Using the AMR code AstroBEAR, we 
have carried out high resolution 3D simulations that follow a stellar mass secondary 
p/| in the co-rotating frame as it orbits a wind producing AGB primary. We first derive a 

ryy resolution criteria, based on considerations of Bondi-Hoyle flows, that must be met in 

order to properly resolve the formation of accretion disks around the secondary. We 
then compare simulations of binaries with three different orbital radii (i? = 10, 15, 
20 AU). Disks are formed in all three cases, however the size of the disk and, most 
importantly, its accretion rate decreases with orbital radii. In addition, the shape of 
^—j | the orbital motions of material within the disk becomes increasingly elliptical with 

increasing binary separation. The flow is mildly unsteady with "fluttering" around 
the bow shock observed. The disks are generally well aligned with the orbital plane 
after a few binary orbits. We do not observe the presence of any large scale, violent 
instabilities (such as the flip-flop mode). For the first time, moreover, it is observed 
that the wind component that is accreted towards the secondary has a vortex tube-like 
fS| , structure, rather than a column- like one as it was previously thought. In the context 

■ of AGB binary systems that might be precursors to Pre-Planetary and Planetary 

^sD 1 Nebula, we find that the wind accretion rates at the chosen orbital separations are 

generally too small to produce the most powerful outflows observed in these systems 
if the companions are main sequence stars but marginally capable if the companions 
are white dwarfs. It is likely that many of the more powerful PPN and PN involve 
s ! ■ closer binaries than the ones considered here. The results also demonstrate principles 

of broad relevance to all wind-capture binary systems. 
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1 INTRODUCTION 

The formation of accretion disks in binary stars can occur 
in several ways. If the orbital separation of the stars is small 
enough then the donor star can overflow its Roche lobe and 
transfer material to the accreting star where conservation of 
angular momentum will force the development of an accre- 
tion disk. This process has been well studied both analyt- 
ically and via simulations (e.g. Eggleton 1983; D'Souza et 
al. 2006). The second mechanism, wind capture via Bondi- 
Hoyle-Lyttleton (BHL) flows, is less well studied and fun- 
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damental questions as to its efficacy remain. In particular it 
is not clear how large an orbital separation is possible that 
still allows for wind capture accretion disks to form. This 
question applies for a range of binary accretion systems but 
our particular focus is on the late stages of evolution for low- 
and intermediate- mass stars. 

Over the last decade the traditional paradigm for plan- 
etary nebula formation, in which successive winds from a 
single star create the nebula, has been challenged on a num- 
ber of theoretical and observation fronts (Bujarrabal et al. 
2001; Soker & Rappaport 2000, 2001; Balick & Frank 2002; 
Nordhaus et al. 2007; De Marco 2009; Witt et al. 2009; 
Huarte-Espinosa et al. 2012; Soker & Kashi 2012). If in- 
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deed binaries play a central role in the formation of pre- 
planetary and planetary nebula (PPN/PN) then we must 
revise our understanding of the evolutionary paths for sun- 
like stars. The basic morphological influence of a secondary 
on the wind poses one set of questions (e.g. Edgar et al. 
2008; Huggins et al. 2009; Kim & Taam 2012ab), but the 
formation of accretion disks around the secondary is of 
more direct importance for constraining which subset of spe- 
cific binary accretion scenarios (Reyes-Ruiz & Lopez 1999; 
Soker & Rappaport 2000,2001; Blackman et al. 2001) can 
lead to the powerful collimated outflows observed in PN 
and PPN (Bujarrabal et al. 2001). Since accretion disks 
plausibly drive collimated outflows via magneto-rotational 
processes jBlandford fc Pavnelll982l;IOuved fc Pudritz |l997l ; 
iBlackman et al.ll200ll ; lMohamed fc Podsiadlowskill2007l ) de- 
termining the conditions for disk formation remains a basic 
issue for studies of the low/intermediate mass star late evo- 
lutionary stages. In particular, understanding the range of 
orbital separations in which disks can form with sufficient 
accretion rates is important for population synthesis models 
that can determine which systems will become PN and/or 
PPN. 

The subject is analytically and computational subtle 
and there have been only a handful of studies addressing the 
aforementioned questions. Thuens & Jorisen (1993) carried 
forward early work on the subject using an Smooth Par- 
ticle Hydrodynamics (SPH) code. Mastrodemos & Morris 
(1998,99; which we will refer to as MM98 and MM99, re- 
spectively) also used SPH simulations to demonstrate that 
disks can form in AGB specific binary systems. Even though 
the resolution was low by contemporary standards they were 
able to set some limits on where disks could form and gain 
insight into a limited set of disk characteristics. In Soker 
& Rappaport (2000) analytic estimates were given for when 
disks would form within AGB binary systems. General stud- 
ies of wind accretion in binaries using a grid based code were 
carried out by Nagae et al. (2004) and Jahanara et al. (2005) 
who explored mass accretion rates and angular momentum 
losses. Recently the issue of disks in AGB binary systems 
has been taken on numerically by de Val-Boro et al. (2009) 
who carried forward high resolution simulations using the 
AMR code FLASH. Their work was carried out in 2-D and 
showed that, in this restricted geometry, disks might form 
out to 40 AU. Mohamed & Podsiadlowski (2007) also stud- 
ied the problem, but using the code GADGET, with a focus 
on the limits of Roche lobe overflow. Disk formation and 
structure have also been explored via semi-analytic models 
by Perets & Kenyon (2012) who found disks forming out to 
100 AU. 

Here we carry out the first fully 3-D high resolution sim- 
ulation study BHL flows and disk formation in AGB binary 
systems since MM98. Our goals are to explore the mecha- 
nisms of disk formation, to study the flow features associated 
with fully formed disks in binaries, and to obtain a first cut 
at the properties of disks formed in these systems. In ad- 
dressing all of these issues, we are specifically interested in 
systems relevant to PPN and PN but our results will be 
broadly applicable to wind capture binaries. 

The structure of our paper is as follows. In section [2l 
we briefly discuss BHL flows in binaries in order to derive a 
resolution condition that should be obtained to accurately 
track the formation of disks in wind capture systems. In 



section [3] we discuss the methods and models used in our 
study. Section [4] provides a description and analyses of the 
results. We then discuss accretion and its implications for 
PPN and PN in section[5] Finally, in section|6]we provide our 
conclusions and their relevance to PPN/PN, BHL studies 
and binary systems as a whole. 



2 RESOLUTION CRITERION FOR BINARY 
BHL FLOWS 

We begin with an evaluation of BHL accretion in the context 
of binary stars. A excel l ent re view of BHL physics in general 
can be found in lEdeari (|2004T ). Our goal in this section is to 
characterize a resolution criterion for tracking the formation 
and evolution of disks in BHL binary systems. 

The general BHL problem occurs when a compact ob- 
ject of mass, m, moves at a constant supersonic relative 
velocity, v r , through an infinite gas cloud which is homo- 
geneous at infinity. The gravitational field of m focuses the 
cloud material located within the BHL radius 



rb — 2Gm/v r 



(1) 



The focused material does not accrete directly onto the ob- 
ject but forms a downstream wake. A conical shock is formed 
along the so called "accretion line" connecting the wake and 
the object. This flow becomes divided into (i) material ac- 
creted onto the object and (ii) material that flows away and 
escapes from the object. These flow components are sepa- 
rated by a stagnation point (Bondi & Hoyle, 1944). 

Here, we study a system in which the cloud material is 
that of an AGB star's wind moving radially outward with a 
terminal speed, v m . The accretor is an orbiting companion 
of mass, m s , with an orbital separation of R > 10 AU and 
orbital velocity, v s . The relative wind speed seen by the com- 
panion is therefore v r = \fv%, + «|. One essential difference 
between the classic BHL analysis and the binary case is the 
accretion flow in the wake. As a consequence of the systems 
orbital motion material, gas at the stagnation point does not 
flow directly toward the secondary but is aimed slightly off 
axis. This occurs because the instantaneous position of the 
companion relative to the BHL accretion flow is constantly 
changing. Hence material in the wake flows towards a point 
between the companion's earlier and current positions; i.e. 
the focused material is accelerated towards a retarded posi- 
tion. Figure [1] shows the considered flow structure. It is, in 
fact, this off-center streaming of the accretion flow (and its 
associated conservation of angular momentum) that allows 
an accretion disk to form. Resolving this process is therefore 
critical to capturing the formation of an accretion disk. We 
note that from the perspective of a co-rotating frame, the 
shifting of the accretion stream occurs due to non-inertial 
forces. 

The time scale, t c , associated with the wind capture 
process scales approximately with the time scale for mate- 
rial to flow through the BHL radius, that is t c — — = 2G 'f 3 , 
where (again) v r is the velocity of the gas relative to the com- 
panion. In an inertial frame moving with the secondary's 
velocity at to, the flow resembles a classical BHL pattern, 
although the secondary is accelerated from rest towards the 

v 2 

primary at o c = The distance, d, travelled by the sec- 
ondary during the wind capture time scale t c is l/2a c t^, 



© 2012 RAS, MNRAS 00Q.ri]fl2l 



Wind- Capture Disks In Binary Systems 3 




Stagnation Region 



Figure 1. Schematic showing the flow structure of the binary 
wind capture process. 
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Figure 2. Plot showing impact parameter as a function of sepa- 
ration for a 1.5Mq primary and 1 Mq secondary. 



tion disk forming or the formation of disk with unphysical 
dynamics such as extreme warps or tilts. 

We note also that our analysis of the impact parame- 
ter implies that disks around secondaries with larger orbital 
separations may be smaller if the initial impact parameter 
is a gauge of the disk final steady state configurations. 



3 COMPUTATIONAL METHODOLOGY 

To model the formation and evolution of disks in binary 
system, we use the adaptive- mesh-refinement (AMR), mag- 
netohydrodynamics code AstroBEAR2.$\. AstroBEAR is 
a highly-parallel, second-order accurate, shock-cap t uring , 
multi-physics Eulerian code (|Cunningham et all 120091 : 
ICarroll-Nellenback et alll201ll ). 
We use BlueHiv<B 

-an IBM massively parallel process- 
ing supercomputer of the Center for Integrated Research 
Computing of the University of Rochester- and 7?anger0 -a 
Sun Constellation Linux Cluster which is part of the Ter- 
aGrid project- to run simulations for an average running 
time of about 0.5 days/orbit using 64-1024 processors. 

We perform calculations in three-dimensions and solve 
the equations of hydrodynamics in a corotating frame, which 
in conservative form are: 

|f + V-(pv) =0 

%i + V • (pvv) = -Vp - pV$ -2p!ixv (4) 
— pfl x (Q x r) 

where p, p and v are the gas density, pressure and flow veloc- 
ity, respectively. We employ an isothermal ideal gas equation 
of state with an adiabatic index of 7 = 1. 

We assume the gravitational potential is due entirely 
to the secondary and has the functional form $ = Gm s /r 
for r ^ 4-^/ dx 2 + dy 2 + dz 2 , and we use a spline softening 
function at smaller radii. While AstroBEAR2.0 has options 
for MHD, self-gravity and heat conduction, we neglect these 
processes in our current study. 



and the distance perpendicular to the accretion line, b, is 
dsina = d^f-. We call b the impact parameter which can be 
expanded as 



b = 



1< / 2Gm s 



\ 2 Vs 

) v r 



2G 2 m 2 3 v 3 s 



(2) 



Solving for v s and r s in terms of the mass of the primary, 
m p , the total mass, M — m s +m p , and the total separation, 
R, we get 



6 = 2 



m?M 3 



U) 



5/2 



1 7/2 



(3) 



where r w = jp$ which for our choice of wind velocity and 
masses gives the profile in Figure [2] 

One of the primary results of the simulations that we 
present here is showing that accurately resolving this phys- 
ical scale in the vicinity of the secondary is essential to cap- 
turing the dynamics of the flow. Since the material initially 
goes into orbit around the secondary with a radius of b, a 
failure to resolve the scale can lead to the lack of an accre- 



3.1 Initial Conditions and Set-up 

Our computational domain consists of a rectangular volume 
of 10x10x5 AU 3 . We employ a fixed grid of 64x64x32 base 
cells plus three, or four (see Table 1), AMR levels of factor 
two refinement. Nested refined cell blocks are centered at 
the secondary's position such that the finest level has 32-64 
finest cells (on the orbital plane) enabling us to resolve the 
central part of the disks and b. 

We employ either wind boundary conditions (sec- 
tion [3TL3]) in the —x,+y and ±z domain faces, or outflow 
only conditions in all other faces (Figure [3]). The positions 
of the primary, the center of mass and the secondary are 
(— r p ,0,0), (0,0,0) and (r 3 ,0,0), respectively. We include Cori- 
olis and centrifugal terms and use a reference frame which 
co-rotates with the secondary. 



1 https:/ /clover. pas. rochester.edu/trac/astrobear/wiki 

2 https://www.rochester.edu/its/web/wiki/crc/index.php/ Sys- 
tems#Blue_Gene.2FP 

3 https:/ /www. xsede.org/web/guest/tacc-ranger 
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Figure 3. Top: Diagram showing the wind model and the grid. 
Bottom: Grid, its coordinates, in AU, and the central mesh. The 
secondary's particle is located at the origin. 

3.1.1 Primary 

We use a simple AGB star model for the primary consisting 
of a mass of 1.5 Mq and a wind with an isotropic velocity of 
v w =10kms -1 and mass-loss rate of m p — 10 -5 Mq yr _1 . 
These values are typical for AGB stars (e.g. Hrivnak et al., 
1989; Bujarrabal et al, 2001). 

We note a few aspects of this model. First, winds from 
AGB stars are driven by radiation pressure acting on dust 
grains. Such interactions are co mplicated and many de tails 
depend subtle non-linear effects l|Sandin fc Hofneril2003l . and 
references therein). Additionally, observations suggest that 
AGB-star winds possess acceleration zones in which the 
rad ial wind veloc ity climbs slowly to its final value (e.g. 
sec 1 ISandin|[200l ). Such acceleration regions may play a 
role in the formation of disks in binaries with separations 
<10AU. We expect this effect to be mild for the sepa- 
rations explored here (10-20 AU), thus no acceleration re- 
gion has been considered and we leave more realistic AGB 
star wind models for future work. Additionally, we note 
that the binary separations that we explore here are out- 
side the wind-Roche-lobe-overflow (WR OLF) limit l|Sandinl 
l2008l : lMohamed fc Podsiadlowskill20ilh . 

3.1.2 Companion 

We model the secondary star using a sink par ticle, based 
on the implementation of lKrumholz et all (|2004 ) . which ac- 
cretes gas within a radius of 4 grid cells. The mass, mo- 
mentum and energy of the gas within this region are then 



acquired by the particle in a conservative fashion. We set 
the particle with an initial mass of 1 Mq in a circular orbit 
about the primary. While important at these radii, we do 
not adjust the orbit for the effect of mass loss or tidal in- 
teractions as they occur on longer evolutionary timescales 
(Nordhaus et al. 2010; Nordhaus & Spiegel 2012). 



3.1.3 Wind solution and injection 

We set the initial conditions throughout the grid with a 
constant temperature of 1000 K and a density given by 

Pw = -r—, t- , (5) 

47r(x p — x) 2 v w 

where rh p —5x 10 _5 MQyr _1 , v w =10kms _1 , x p and x 
are the positions of primary's orbit and that of an arbitrary 
grid cell relative to the center of mass (0,0,0), respectively. 
We calculate the velocity field of the solution by solving for 
the characteristics that leave the surface of the primary at a 
retarded time, t r = t—\x.\/v w , with a velocity vector pointing 
towards x. Wind speeds are chosen so that v w > |v s |, where 
v s is the secondary's orbital velocity. We assume that the 
distance from the primary's surface to x is larger than |x p |, 
which constrains the separation between the secondary and 
the grid's boundaries. As time goes from t r to t, the primary 
covers a circular segment of radius r p starting at x p (i r ), with 
a displacement vector d = x — x p (i r ). We calculate the wind 
normal, n, so that (v m ii + v p (t r )) x d = 0. Starting at the 
primary's surface, the wind velocity field, u = u(x,t), is 
then 

u = v w h + v p (i r ), (6) 
which yields a better approximation of the retarded time 

r r = t-|x|/|u|. (7) 

We iterate the above computations until we find a con- 
vergent wind solution for the gird cell located at x. Finally, 
we transform the solution to a reference frame which co- 
rotates with the secondary. 

Once the initial conditions are set, we continually give 
the above wind solution to the grid cells in the — x, +y 
and ±z domain faces. For each iteration however we use 
t t , x p (r r ) and u as initial conditions in a recursive Runge- 
Kutta method to calculate the deflection of u -as it travels 
from the primary's surface [~ x p (r r )] to x- caused by the 
secondary's gravity field. We consider the orbital motion of 
the secondary in these computations. This step yields yet 
better estimates of both n and r r than the calculations in 
equation [6] alone. We therefore use the new values of n and 
r r for the next iteration of the wind solution. 

Finally, to match the initial conditions and the injected 
wind solution we allow the gravitational effect that the sec- 
ondary has on the gas to increase linearly during one wind 
crossing time, l.6rb/v w , from zero to Gm s /r. The field re- 
mains constant thereafter. This has no effect on the binary 
orbital motion. 
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Table 1. Simulations and parameters. 



Name Separation [AU] Resolution r b [AU] b [AU] 

Ml 10 64x64x32+3amr 4.9 2.4 

M2 15 64x64x32+4amr 5.8 1.5 

M3 20 64x64x32+4amr 6.3 1.0 



tours of the disk in slices though the orbital and perpendicu- 
lar planes. Each panel has 12 contours arranged by line type 
and color: blue, red and black represent densities of 3.5, 
10.5 and 35 xlO 10 partem -3 , respectively. Dashed, dotted 
and dashed-dotted lines represent densities at t = 1, 2 and 
3 orbits, respectively. 

The density contour maps for different orbital separa- 
tions are mutually consistent in terms of evolutionary stages 
but are correspondingly different with respect to the time 
scales and some features of these phases. Comparing con- 
tours of the same color (Figure [6| we see a given iso-density 
structure at different evolutionary times. These show that 
once the disk forms further structural changes are mild. 
Thus in all cases the disks begin (t — 0.1) from long elliptical 
structures that can be identified as the curved BHL accre- 
tion stream. By (t = 1.0) the inner regions of the disk have 
formed and we see well defined inner black contours which 
do not change appreciably with time. The middle and outer 
contours in all three case (red and blue) do show some evo- 
lution after t = 1.0, however these changes tend to be higher 
order in morphology whereas overall shape (circular in the 
a — 10 and 15 AU case and elliptical in the a — 20 AU case) 
and radius do not change. 

Comparing the solid lies in left panels we find, in agree- 
ment with the discussion in sectionfJl that the angle between 
the accretion flow and the wind decreases with a. This angle 
is close to zero already in the 20 AU model. Thus, as ex- 
pected, the solution ([2| converges to the BHL one for large 
separations. 

We note that given our choice of 7 = 1, the accretion 
disks which form in our models are thin which is consistent 
with a low gas temperature and sound speed (h/R ~ c/Sl). 
The disks orbits are also not circular. We see orbits with 
eccentricities e > in all cases however the magnitude of 
average eccentricity increases with increasing orbital separa- 
tion a. And while our disks are thin, we do see flared vertical 
structures that are also not symmetric in the (r, z) plane of 
the disk but larger on their downstream side relative to the 
wind. The effect of the wind on the disk is also apparent 
in that the steepest density gradients are located where the 
disks impinge on the incoming stellar flow (top left region 
in left panels and left region in right panels, Figure 0). We 
find that both the disks' outer radius, r^, and edge height 
are inversely proportional to a. In section |4~31 we show that 
the shape of the gas parcel orbits in the disk is a function of 
a as well. We note that the disk inclination angles are gen- 
erally small meaning the disks are close to perpendicular to 
the orbital plane. However small time variation in the flow 
occur even after the disks settle into a quasi-steady state. 

4.2 Disk mass 

In Figure [7| we show profiles of the disks' mass, Md(t) (solid 
lines), as a function of time and binary separation. These 
are calculated by summing the mass of grid cells containing 
bound gas (those having gravitational energy greater than 
kinetic energy). 

In the 10 AU case we see that the disk mass in- 
creases during the first orbit and then saturates at about 
7xlO" 6 M Q . The disk mass then oscillates with a period of 
~1 orbit and an amplitude of about lxlO~ 6 M0. 

To understand the origin of the disk mass oscillations, 



4 RESULTS 

We carry out three simulations corresponding to stellar sep- 
aration of 10, 15 and 20 AU. Table 1 summarizes the models 
and their relevant parameters. 

4.1 Disk structure and formation 

In order to present the main results of our study in Fig- 
ure [4] we show logarithmic density gray-scale and velocity 
field maps of our three simulations at t — 3 orbits. Rows 1 
and 2 of Figuref4]show slices though the orbital plane (at two 
different scales) , while the bottom row shows slices though a 
longitudinal plane that intersects the secondary. These den- 
sity maps allow us to see the large scale structure of the 
flows and to compare those flows as the orbital separation a 
is increased. 

First we note the basic BHL converging flow which can 
be seen in bottom row of Figure 2] The AGB wind can be 
seen streaming in from the left and is subsequently focused 
towards the mid plane by the gravitational force of the sec- 
ondary. At the time of these images the global flow is con- 
siderably more complicated than that described by the sim- 
ple BHL scenario however because of the presence of fully 
formed accretion disks. 

Thus, as Figure [3] makes clear, accretion disks are cre- 
ated in all our simulations. These disks form, as expected, 
from the BHL accretion stream that originates at the stag- 
nation region downstream from the secondary. We provide a 
better illustration of the stagnation region in Figure [5] where 
the top left panel shows iso-surfaces of the disk's density and 
slow material's velocity. The top right panel shows one flow 
streamline which starts at the wind, is deflected by the sec- 
ondary, reaches the stagnation region then, and is finally 
accreted onto the disk. In the bottom panel we show sev- 
eral flow streamlines from a different perspectives using a 
velocity color table that clearly shows the stagnation region 
structure. We note the appearance of nested flow surfaces, 
whose areas increase with distance from the stagnation re- 
gion. Thus the accretion column resembles a strong vortex 
tube and our simulations are the first to capture this aspect 
of BHL flow in binaries. 

The density maps in both the equatorial and perpen- 
dicular planes (fig [4]) also reveal the presence of bow shocks 
around the accretion disks. Such shocks also form in the 
classic BHL flow however once the disks emerge, incident 
AGB wind must stream around them creating more complex 
post-shock flow patterns. The velocity vectors in the upper 
row panels of Figuref3]are of particular interest, showing the 
redirection of high density flow behind the bow shock and its 
interaction with the counterstreaming flow in the accretion 
disk. 

To understand the structure and evolution of the disks 
more explicitly in Figure[6]we show logarithmic density con- 
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Figure 4. Logarithmic density gray-scale and velocity field maps for the 10 AU (left), 15 AU (middle) and 20 AU (right) models at 
t = 3 orbits. The primary is simulated to be to the left of the maps. The particle is receding in bottom row panels due to its orbital 
motion. 



we calculate the gas flux leaving the grid though both the +x 
and — y boundaries (see Figure [3]) in Mq yr _1 (dashed lines, 
Figure 0. The position and shape of the gradients in M<j(£) 
compared with the mass fluxes indicate that the disk's mass 
decreases with increasing outbound flux and vice versa. This 
suggests the wind's post bow shock ram pressure mediates 
the disk mass; the wind appears able to strip a small fraction 
of disk's gas and liberate it from the gravitational pull of the 
companion. 

Moving to comparisons with simulations with different 



a, we see that the disk mass in the 15 AU model reaches its 
maximum value of about 4x 10 -6 Mq in ~ 0.7 orbits (ear- 
lier than the a — 10 AU case). Once again we see disk mass 
oscillations with amplitudes of about lx 10~ 6 M© and a fre- 
quency of ~ 0.3 orbits. Note that these oscillations have a 
smaller amplitude and period than for the 10 AU case. This 
would be consistent with our wind stripping interpretation 
as the wind's ram pressure is proportional to a~ 2 x^ 1 , where 
xj 1 is the distance between the secondary and the center of 
mass. 
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Figure 5. Structure of the stagnation region for the 15 AU case 
at t = 3 orbits. Top: gray disk density iso-surfaces and blue slow 
material velocity iso-surfaces. Middle: flow streamline that starts 
at the wind, is deflected by the secondary, reaches the stagna- 
tion region and is finally accreted onto the disk. Bottom: flow 
streamlines in Mach units. 

The disk mass profile of the 20 AU model is quite differ- 
ent. While an initial disk forms with mass ~6.5xlO~ 7 M Q 
a sharp growth in mass is seen at t « 0.85 orbits. The in- 
crease occurs in less 0.1 of an orbit (t — 5.5 yr) and the disk 
mass then reaches a maximum value of ~ 6.5 x 10" 7 M af- 
ter a brief period of relaxation. Note that in this case we 



do not see any disk mass oscillations. This behavior is con- 
stituent with our wind stripping interpretation as the wind 
ram pressure is the weakest in the 20 AU case. 



4.3 Disk orbits 

An interesting result from our simulations is that the shape 
of the disk gas orbits significantly depend on a. We show 
this relation using the disk gas streamline maps in Figure |8j 
where the top panel presents a 3-D map of the disk gas 
velocity distribution in colors for the 10 AU model, along 
with sliced (half the disk) density contours in translucent 
red color. The lines show that the gas is supersonic, the 
orbits have low eccentricities, e, and the velocity decreases 
with distance from the particle. 

In the bottom row of Figure [8] we show a quantita- 
tive comparison of disk gas orbits as a function of a, at 
t =3 orbits, corresponding to the 10, 15 and 20 AU cases 
from left to right, respectively. The maps consistently show 
that e is a function of both, and a. For a particular disk, 
material located at small disk radii shows moderately higher 
e than material located at larger radii. Comparing disks 
with different a we see that (i) the magnitude of average 
e increases with increasing a and (ii) the angle between the 
semi-major axis of the orbits and the wind velocity field de- 
creases with a. The latter is particularly clear in orbits with 
small disk radii. This behaviour is in good agreement with 
our findings in section HJ] the angle is close to zero in the 
20 AU model (right-most panel). 

Finally, the color bars in the bottom panels of Figure [5] 
show disk gas position along the direction perpendicular to 
the orbital plane; receding material has red color whereas 
approaching material has blue color. Note that we have cho- 
sen different color limits in each panel in order to stress disk 
orbital features. The color of the lines in the 10 and 15 AU 
cases (left and middle panels) consistently shows that mate- 
rial oscillates about the orbital plane mildly with amplitudes 
no larger than 0.2 disk radii and frequencies of order half an 
orbit. The oscillation are much weaker in the 20 AU case. 
We do not see any correlation between the color of the lines 
-i.e. the position about the orbital plane- and the wind ve- 
locity field direction (which points from the top left corner 
to the bottom right one; Figure [8] bottom) . 



5 DISCUSSION 

5.1 Accretion onto the companion and 

implications for PPN/PN phenomenology 

In Figure [9] we show the evolution of the sink particle's 
accretion rate. The plots show the total gas flow into a region 
within a 4 cell radius from the particle's center (section l3.1.2|l 
as a function of time. The plots demonstrate that after a 
steep initial increase which lasts for about 0.05 orbits the 
accretion rate reaches a quasi-steady value with oscillations 
of a period that decreases wi th a. These accretion histo ries 
are consistent with those of Ide Val-Borro et al.1 (2009, sec 
their Figure 12). We note that we see a small dip in the 
20 AU profile, close to 0.9 orbits, which corresponds to the 
disk formation time. 



© 2012 RAS, MNRAS 000,rflfl2l 



8 M. Huarte-Espinosa et al. 



1 orb»J(rflci, 1 orb-dgsh, 2orb"dot, 3 orb-dqih-dot 




3.6; 



2.0 



1.0 



0,0 



.1.0 



-2.0 



-3.0 



4b 






3« 10" 16 


"(1,3. ICS 












/ 




T — ' 


> 










1. 1 

V r l' 
i 

H- 


<■- 

! 


4 




hi 
\ 








* 4. 


1 
h 




- 1 

b 










% 




H i 


r 







































-3.0 -2.0 -1.0 0.0 1.0 2.0 3.0 



Soiti-irjhd, laib»dniri, SoiC-dnr, inib'Snih iJc' 
<blu». rod. bract) - i.J«1D« Id ■ (J. 4 IMS pan/cc 





' , -lo-"'-asr^ra =s= ffr = T 



1.0 
0,3 
0.0 

-0.5 
-1.0 



■ 




































• 1 IT 1 . „ 


















-1.0 -1 


.0 o 


1 


o a 


o i 






Air 




2.0 



Figure 6. Disk density structure at 4 times. Left: face-on view where the wind enters the grid from the top left corner at an angle close 
to 45° and the particle orbital motion is upwards. Right: edge-on view where the wind enters the grid horizontally from the left side and 
the particle is receding due to its orbital motion. The primary is simulated to be to the left of all maps. Top, middle and bottom rows 
correspond to binary separations of 10, 15 and 20 AU, respectively. 
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The analytic Bondi-Hoyle accretion rate for a high ve- 
locity wind and negligible pressure conditions is given by 



(8) 



where Rt =2Gm s (y^ + Vg) 1 is the modified BHL ra- 
dius (Edgar, 2004). For our 10, 15 and 20 AU models, 
m B H take the values of 3.3, 1.8 and 1.2x 10~ 6 M© yr _1 , 
respectively. Thus the secondary accretion rates we find 
(Figure [9]) are m 10 Au ~0.9m S H,io, rhisAU ~0.2m S H,i5, 
rri20Au ~ 0.16msH,2o, where the subscripts BH, 10, BH, 15 
and BH, 20 indicate the Bondi accretion rate at 10, 15 and 
20 AU. 

Our results appear to consistent with those found 
by other author s . Fo r example in the models of 
Ide Val-Borro et"al] (120091 ). a of 70 AU separation produced 
Macc/M w i n d ~0.06 after 2 orbits. While we have not run 
simulations with such a large orbital separation, these val- 
ues are in line with the trend we observe in our results. We 
note that in general we should somewhat expect higher ac- 
cretion rates in our 3D simulations than in their 2D models, 
since only 3D runs can capture poloidal accr etion. The semi- 
analytic models of iPerets fc Kenvonl (|2012l ) yield m values 
that are also consistent with ours. For example their Model 
2 with a = 10 AU, m p = 3.0 M ot j t, m s =1.0M o( j o t and 
M p = 10" 5 Modoty" 1 yielded M„ cc ~ 3xl0~ 7 . We note that 
the models of lPerets fc Kenvonl (|2012l ) are semi-analytic and 



so a direct comparison between our work and theirs is not 
straightforward . 

The results shown in Figure [9] are relevant for constrain- 
ing the class of models that can power PN outflows and, 
in particular, the PPN outflows in the reflection nebulae 
thought to precede the PN phase. The observations of Bu- 
jarrabal et al. (2001) imply that statistically, many PPN 
seem to harbor very powerful outflows whose energies and 
estimated acceleration times give mechanical luminosities as 
high as 4 x 10 36 ergs _1 (Blackman 2009, see also Huggins 
2012). Can the accretion rates we find power such high lu- 
minosities? 

If these high powered PPN outflows do indeed result 
from accretion disk around a secondary, then one typically 
expects no more than 10% of the accretion rate to turn into 
mass outflow and the outflow speed fir a = Qvk (e.g. Frank 
& Blackman 2004) where Q is the orbital speed the inner 
edge of the disk, ta is the radius of the Alfven surface, vk 
is the Keplerian speed at that radius and Q is a numerical 
factor determining the jet speed as a multiple of the Keple- 
rian speed (if the jet propagated into a vacuum). Although 
PPN jets would propagate into the stellar envelope such that 
the observed jet speed would be reduced using momentum 
conservation (Blackman 2009), the mechanical luminosity 
of the jet should be comparable to its naked value if energy 
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Figure 8. Top: 3-D disk gas orbit streamlines and sliced density contours (red) for the 10 AU case at t =2 orbits. The central white 
small sphere is the secondary's particle. Streamline colors denote orbital speed in Mach units. The wind enters the grid towards the 
image. Bottom: Disk gas orbit streamlines at t =3 orbits corresponding to the 10, 15 and 20 AU cases, from left to right. Grid squares 
arc 1 AU 2 . The wind enters the grid from the top left corner moving towards the bottom right. Colors denote distance perpendicular to 
the orbital plane in AU and indicate the orbits inclination. 



is conserved in the outflow. Accretion would provide a jet 
mechanical luminosity of 



1 Mac n2„,2 
2~L0~ { * V K 



= 4 x i(r 



5. 6xl0- 5 M /yv 



(f) 2 ( 



V K 



440km/ b 



(9) 



where we have scaled vk to that of an orbit at the surface 
of a main sequence star of 1 Mq and 1 R@ and taken an op- 
timistic Q = 3. The accretion rates we find (see Figure [§]) 
would then fall short by ~ 2 orders of magnitude if the sec- 
ondary is a main sequence star. However, for a white dwarf, 
v% is ~ 70 times larger and the required accretion rate to 
produce the same mechanical luminosity is reduced by the 
same factor, and becomes 8 x 10 -7 Mq yr _1 . This rate is 
not that dissimilar from the accretion rate we found for the 
10 AU separation case. Note that the estimates above de- 
pend on the value of Q which is somewhat uncertain and 
on the actual inner most radius of the disk used to compute 
vk- It could be that the disk does not extend all the way 
to the star in which case vk in the above formula would be 
reduced. 

The fact that only a white dwarf companion could sup- 
ply the outflow powers for the highest powered PPN, high- 
lights the need to determine the mechanical luminosity dis- 
tribution function of PPN. If the highest outflow powers 
are much more common than the expected frequency of 
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Figure 9. Secondary (particle) accretion rate evolution. 



white dwarf companions within 10 AU then the required 
mechanism would likely require closer binary interactions 
and perhaps even common envelope evolution. More work is 
needed to assess the various modes of accretion that can 
operate within a common e nvelope, either onto the pri- 
mary, or onto the se condary l|Nordhaus fc Blackmanll2006l : 
iNordhaus et al.ll201ll ). It should be noted that the late stage 
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PN are less demanding with respect to required accretion 
rates than the early PPN objects and it may be that not 
all PPN are as demanding as the Bujarrabal et al. (2001) 
sample. 

One example of a lower power object of interest is the 
Red Rectangle PPN. In fact this object, reported in Bujarra- 
bal et al. (2001) is one of the lower powered PPN presented 
therein. Witt et al. (2009) showed that the Red Rectangle 
is consistent with accretion onto a main sequence compan- 
ion at binary separation <1 AU. The jet has a mechanical 
luminosity ~ 7x 10 33 ergs~\ Witt et al. (2009) plausibly 
argue from presumed evidence for disk emission that the disk 
around the companion is accreting via Roche lobe overflow 
from the primary at a rate 2-5x 10 -5 Mq yr _1 . However, 
these authors argue that the jet outflow speed would be 
launched at 100 km s" 1 given this accretion rate, which is 
significantly below the escape speed at the inner edge of the 
disk if the disk were to extend all the way to the star. Such 
an accretion rate could certainly power the observed jets in 
the Red Rectangle. 

In the future, a larger parameter survey covering a sys- 
tematic range of grid resolution setups and accreting radii 
about the sink particle -gravity softening radii- would be 
helpful both, to further constrain the connection between 
core accretion and jet launch in binary-formed disks and to 
better characterize numerical artifacts. 



6 CONCLUSIONS AND SUMMARY 

We have performed the highest resolution simulations to 
date of wind accretion and disk formation in an evolved star 
binary system. We simulated the wind accretion process at 
three different orbital radii, 10, 15 and 20 AU, and found ac- 
cretion disk formation in all cases. We have also derived an 
impact parameter criterion that determines the minimum 
required resolution for such simulations to produce phys- 
ically consistent results. This impact parameter measures 
the perpendicular (retarded) distance from the stagnation 
point in the BH accretion stream to the accretor and as 
such is a measure of the distance at which material from the 
wind is captured into orbit around the secondary. We find 
the impact parameter decreases non-linearly with increasing 
orbital separation. 

Additionally, accretion disks formed in our at 10 
and 15 AU simulations incur mass oscillations which appear 
to result from ram pressure stripping of the outer disk radii. 
The 15 AU oscillations become weaker as the wind ram pres- 
sure decreases with orbital radius. The 20 AU disks, showed 
no oscillations, consistent with an even weaker wind ram 
pressure. 

The disk accretion rates evolve to be quasi-steady af- 
ter 1/4 orbits and decrease as a function of orbital radius, 
ranging from ~3xlO -7 M yr~\ ~ 5x 10~ 7 M yr" 1 and 
SxlO^Moyr -1 for the 10, 15 and 20 AU cases, respec- 
tively. The standard Bondi-Hoyle model of accretion applies 
only at large separations in that only at such separations 
is the angle between the accretion flow and the wind zero. 
We see this angle becoming < 5° for the 20 AU case but it 
is non-negligible for the smaller radii runs. While the basic 
problem of accretion from a wind studied herein has broad 
applicability in a number of contexts, we chose the parame- 



ters to be consistent with those of a stellar mass companion 
accreting from an AGB wind. This is motivated by the need 
to identify constraints on the mode of accretion that might 
be needed to explain the accretion powered jets and out- 
flows of PPN and PN. Toward this end, we find that the 
accretion rates from our simulations fall short of what is 
needed to be able to power the most powerful young bipo- 
lar PPN unless the companion is a white dwarf. Closer bi- 
nary interactions, which includes common envelope evolu- 
tion, are likely an important player for some PPN objects 
but; our present work focused only on binary radii at 10 AU 
and larger separation in order to be be largely outside of 
the radiative wind acceleration zone. It is useful to compare 
our work with other studies of wind capture accretion in 
binary systems. Theuns & Jorissen (1993) used SPH simu- 
lations and despite low resolution by today's standards were 
able to capture the basic features of a bow-shock and disk 
formation. In the seminal work by Mastrodemos & Mor- 
ris (1999) SPH methods were again used to systematically 
explore the formation of accretion disks. Their work implied 
that disks could form out to distances of at least 21 AU. 
Their 3-D models also hinted at the formation of signifi- 
cant warps in the disks which would be important as these 
might drive p recession in any jets which form from the disks. 
The work of Ide Val-Borro et all (|2009r ) found disks form- 
ing out to radii of 70 AU. If such a result where to hold it 
would be important for establishing the population of bina- 
ries fca_wdiich_cUsl£_wmds_jmght power jets. Unfortunately 
the Ide Val-Borro et "all |2009) study was carried out in 2-D 
only so it not possible with our work to confirm their outer 
limits on disk formation. Of equal importance was the role of 
the inner AGB wind laun ching radius R w in their studies. 
Ide Val-Borro et all (|2009l ) allowed the inner wind launch- 
ing radius to be a free parameter. The justification for this 
was the possibility of large wind acceleration z o nes fo r AGB 
stars. In their simulations Ide Val-Borro et all l|2009h found 
that as R w increased the flow pattern became more complex 
appearing as a mix between Roche Lobe overflow and wind 
ca pture. Such a mix of accret i on m odes was also implied 
by iMohamed fc Podsiadlowskil (|2007h . In iPerets fc Kenvonl 
l|2012l ) semi-analytic methods also found disks forming at 
large binary separations (a = 100 AU). If these results can be 
verified in full 3-D fluid dynamical simulations they would 
also enlarge the population range of outflow bearing sys- 
tems (though these outflows would likely be of low power as 
explained in the previous section). Taken as a whole under- 
standing these possibilities in greater detail will be essential 
for understanding the evolutionary pathways and outflow 
generation of PPN and PN systems. 

Finally we note that our study provides new insights 
into the Bondi-Hoyle accretion in binary systems. We find 
for example that the stagnation region in which the accretion 
"line" separates into flows directed towards and away from 
the accretor is best viewed as a vortex tube as instream- 
ing material flows along helical paths towards the secondary 
star. In addition we find the flows (and most important the 
disk s) which form via binary BH flows appear qui te sta- 
ble (|Blondin fc Popdl2009l ; iBlondin fc Ravmerll2012r i. Only 
small amplitude fluctuations in the orientation of different 
disk annuli are seen and these would appear to be too small 
to drive significant jet precession. In conclusion, more work 
is needed to explore the consequences of interacting binaries 
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for a wider range of orbital radii. And more observations are 
needed to determine the mechanical luminosity function of 
PPN jets to assess what fraction have the high powered jets 
that are typical of the Bujarrabal et al. (2001) PPN sample. 
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